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ABSTRACT 

Type la supernovae (SNe la) have relatively uniform light curves and spectral evolution, which make 
SNe la useful standard candles to determine cosmological parameters. However, the peak brightness 
is not completely uniform, and the origin of the diversity has not been clear. We examine whether 
the rotation of progenitor white dwarfs (WDs) can be the important source of the diversity of the 
brightness of SNe la. We calculate the structure of rotating WDs with an axisymmetric hydrostatic 
code. The diversity of the mass induced by the rotation is ~ 0.08M Q and is not enough to explain the 
diversity of luminosity. However, we found the following relation between the initial mass of the WDs 
and their final state; i.e., a WD of smaller initial mass will rotate more rapidly before the supernova 
explosion than that of larger initial mass. This result might explain the dependence of SNe la on their 
host galaxies. 

Subject headings: binaries: close, supernovae:general, white dwarfs 



1. INTRODUCTION 



Type la supernovae (SNe la) are widely believed to be 
thermonuclear explosions of mass-accreting white dwarfs 
( WDs ) , whose masse s have r eached a certai n critical level 



(e.g.. lNomoto et all l|1997|) : iBranchl (|1998f> for reviews). 
However, the observed brightnesses of SNe la have a 
little diversity. SNe la have been standardized by ap- 
plying the empirical relation existing between the peak 
brightness and the li ght curve width which is obtained 
from nearbv SNe la fPhillipd 119931 lHamuv et aD 1199(1 
IRiess. Press fe Kirshnerlll996j) . Explaining this diversity 
is very important to the use of SNe la as standard candles 
for the determination of the cosmological parameters be- 
cause whether the relation is valid for high-redshift SNe 
la is still uncertain. The diversity is likely to be related 
to the evolution of the progenit or systems l)Branchlll995t 
lLivioll200(l iNomoto et aT1l2f)O0ft . 

Theoretically, the critical mass for the explosion of 
SNe la depends on their progenitor models. There are 
two groups of models according to the masses of ex- 
ploding WDs: Chandrasekhar-mass models and sub- 
Chandrasekhar-mass models. In the Chandrasekhar- 
mass models, a WD is supposed to explode at the cen- 
ter when its mass becomes very close to or exceeds the 
Chandrasekhar's limiting mass (~ 1AM@). Because of 
presence of a limiting mass, the property of the WDs is 
almost identical in the Chandrasekar-mass model. This 
property would lead to the almost uniform peak bright- 
nesses of SNe la. In the sub-Chandrasekhar-mass mod- 
els, on the other hand, an explosion is initiated at the 
bottom of the helium layer before th e WD mass reaches 
the C handrasekhar limiting mass l|Wooslev fe Weaverl 
1994). Inward deflaglation waves converge at the cen- 
ter of the WD and ignite detonation waves, whic h 
blow off the whole WD fWooslev. Taam fc WeaverF l9861. 



Sub-Chandrasekhar-mass models have spectra that are 
too blue and are not in good agre ement with the ob- 
served spectra of s tandard SNe la ( Ho flich et alJ Il99"rl 
iNugent et aHll997h . However, some of the peculiar SNe 
la, su ch as SN 2000 cx have bluer spectra than standard 
ones i|Li et aUl200H) and may be explained by the sub- 
Chandrasekhar-mass models. 

Two possible scenarios have been considered for the 
Chandrasekhar-mass models. One is the double degen- 
erate (DD) scenario in which two WDs in a binary sys- 
tem merge together as gravitational r adiation extracts 
the orbital energy f rom the system l)Iben fc Tutukovl 
119841: lWebbmklll984|) . As a result, the two WDs turns 
into one WD and a thick disk surrounding it whose to- 
tal mass exceeds the Chan drasekhar mass. However, 
iSaio fc Nomotol l)1985lll998j) shows that the merged WDs 
cannot result in SNe la because of the relatively high ac- 
cretion rate from a thick disk. They are likely to turn into 
O+Ne+Mg WDs, then finally become neutron stars via 
accretion-induced collapse. The other is the single degen- 
erate (SD) scenario in which a C+O WD in a close binary 
system accretes hydrogen- rich matter from its companion 
star to grow to the Chandrasekhar mass. Such systems 
are actually observed as supersoft X-ray sources. There- 
fore, Chandrasekar mass model in the SD scenario is the 
most favored model now for standard Sne la. 
The evolution of accreting WDs have been extensively 



investigated by many aut hors iNomoto et alJ 11982a : 
I Kohavashi eHnligM Ettifl ILi fc van den HeuvelfTl997 : 
lLangerll2000() . To make a supernova explosion, a WD 
needs to get enough mass from its companion star. When 
a star becomes a WD in a binary system, its companion 
star is still on the main sequence. As the companion 
evolves, its envelope expands and fills its Roche lobe. 
Then gas flow toward the WD starts and forms an accre- 
tion disk. The WD begins to gain the mass and angular 
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momentum from the rotating disk. 

The evolution of a mass-accreting WD depends mainly 
on the accretion rate (M). A W D needs an appropriat e 
accretion rate to become a SN la l|Nomoto et al.l lT982b). 
Too low M (M < 10 _8 M Q yr _1 ) causes a surface nova 
explosion and the WD loses its mass. Too high M causes 
the WD to expland to the size of a red giant and leads to 
the formation of a common-envelope. After the com- 
mon envelope forms, the gas that has overflowed the 
Roche robe of the system reduces the angular momen- 
tum of the system, which induces merging of the two 
stars l|Nomoto et alJlT97fl|) . It has been thought that a 
SN la is not expected to occur for such very low M or 
very high M. 

However, the strong wind solution is found for M over 
the certain critical value using updated OPAL opacity. 
If M > Af wm j the WD blows stro ng winds to av o id the 
comm o n-envelope forma tion (see lHachisu et alJ l)1996l 
1999a); H achisd <Jl999b) for details of the wind model). 
Therefore, more WDs can accrete mass from the com- 
panion stars until the mass reaches the critical mass for 
SNe la. 

Previous studies on the SNe la progenitors have ig- 
nored rotation or used averaged centrifugal force as - 
suming spherical sym metry (e.g., lLangerl l)2000l |2003|) ; 
iPiersanti et al.1 l)2003ft ). However, it is natural to think 
that accreting WDs are rotating because they must gain 
not only mass but also angular momentum from the ac- 
cretion disk. If rotation of the accreting WD is fast 
enough, the shape of the WD is distorted and the dif- 
ference from the s pherical case cannot be ignored (e.g., 
IWamret all lEffll h 

Although it is very important to consider rotation of 
the progenitors properly, there has been little system- 
atic study on the multi-dimensional structure of rotat- 
ing WDs in connection with the brightness of SNe la. 
Our main objective are to investigate the strucuture and 
evolution of accreting WDs in close binaries with a two- 
dimensional code and to see (1) to what extent the ro- 
tation of the WD affect the strucuture and (2)to what 
extent the diversity of SNe la due to the progenitor's 
rotation can be obtained. 

An accreting WD gains angular momentum from the 
rotating disk as well as mass, and gets rotating faster and 
faster. Rapid rotation affects the mass and structure of a 
WD just before the supernova explosion. The existence 
of centrifugal force effectively reduces gravity and thus it 
is possible to sustain a larger amount of matter with the 
same degenerate pressure. Thus, the critical mass will 
increase as rotation becomes faster. 

When the rotating WD reaches the critical mass (to 
some extent larger than that of nonrotating WDs), it 
will explode as a SN la. The increase of mass and the 
change of inner structure will affect its nucleosynthesis. 
The brightness of a SN la is mainly dependent on the 
amount of 56 Ni synthesized in the explosion; therefore 
rotation of progenitor WDs will affect the brightness, too. 

In §2, we describe models and assumptions. The re- 
sults of our calculations are shown in §3. Disucussion is 
in §4 and the conclusions are summerized in §5. 

2. MODELS 



2.1. Basic Assumptions 

To calculate the structure of rotating WDs, we use 
Hachisu's axisymmetric code (Hachisu Self-Consistent 
Field [HSCF] method: IHachisul twM) ). The number 
of meshes are 129 in both radial and angular directions. 

We adopt the equation of state of complete degenerate 
electron gas for whole WD. 



P = A 



p = Bx 3 



x(2x 2 - 3)(x 2 + I) 1 ' 2 + 31n(V 1 + x 2 + x) (1) 

(2) 



where 



PF 



Pf 



-z—rie , n e = (3) 

87T J \i e 

A = 6.002 x 10 22 dynes • cm" 2 (4) 



B = 9.736 x 10 5 /i e g • cm" 



(5) 



and /i e is the mean molecular weight of electrons and 
fixed as /i e = 2 in all calculations.; No is the Avogadro's 
number, n e is the mean electron number density and 
Pf is the Fermi momentum. No Coulomb interaction is 
included. 

We assume that a WD is at rest when it starts to 
gain mass and angular momentum; i.e. J = Q = at 
M = Mint, where M- mt is the initial mass of the WD. 
Angular momentum of accreting gas is supposed to obey 
Keplerian rotation law, and gives all of its angular mo- 
mentum to the WD. Under these assumptions, specific 
angular momentum j is given as follows: 



3 — jKepler 
^surface = ~G l 



surface? 

P(r') 



■d 3 r- 



(G) 
(7) 



where R is the equatorial radius of a WD and "^surface 
is its surface gravitational potential. At each step of the 
calculation, the WD is kept in hydrostatic equilibrium by 
the balance between the gravitational force, centrifugal 
force, and the pressure gradient. 



P 



1 dP + ^ 



Q 2 wdw = C(constant), (8) 



where f2 is the angular velosity and w is the distance 
from rotation axis; i.e., w = r(l — p 2 ) 1 ^ 2 in spherical 
polar coordinates (r, fi,(f>). 

The mechanism and timescale of the angular momen- 
tum transfer in the star have not been well understood. 
We thus adopt the assumption that the WD is separated 
into two regions. One is the central, nonrotating region, 
and the other is surrounding rigidly rotating region, al- 
though i n reality there exists a gradient of the a ngular 
velocity l|Yoon et al.ll200l lYoon fc Langerll2003lK The 
ratio of the mass of the former region to the whole WD 
is denoted as, 

( = ^, (0) 

where Mnr is the mass of the nonrotating region. A 
WD is initially at rest, i.e., £o = 1- We consider two 
extreme cases of anglar momentum transfer: (1) Angular 
momentum transfer is fast and the whole WD is rotating 
rigidly (£ = 0). (2) Angular momentum transfer is slow 
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and initially nonrotating WD remains static while the 
accreting gas is rigidly rotating around it (£ = Mj nt /M). 

In each case, the angular velocity fl is expressed as 
follows: 

il = n (= constant) (M(w) > M NR ) (10) 
fl = (M(w)<M NR ) (11) 

where M(w) is the mass included in radius w. 

The HSCF method is a variation of the self-consistent 
method to solve the axisymmetric stellar struc ture and 
chara cterized by the selection of parameters llHachisul 
1986). To obtain the structure of a WD, one should 
specify its central density (p c ) and the axis ratio (q) in 
the HSCF method. The axis ratio (q) of the WD is the 
ratio of its polar radius to its equatorial radius, which 
is an index of the strength of rotation. If a WD rotates 
rapidly, its shape is distorted by the centrifugal force and 
its polar radius becomes shorter than its equatorial ra- 
dius. For a spherical WD, q = 1, while, a WD with 
higher angular velocity has smaller q. 

In this calculataion we need another parameter £ (the 
ratio of the mass of the nonrotating region to the whole 
WD; ea 1101 ITT]l to specify the structure. We have calcu- 
lated more than 100,000 WD models with different sets of 
(p c , q, £). For each £, the ranges of p c and q are bounded 
by the criteria of the supernova explosion and the critical 
rotation. Each criterion is described in §2.2. 

2.2. Criteria 

The WD is compressed as its mass increases. On the 
other hand, rotation mitigates compression by reducing 
effective gravity. For given J, the maximum density /3 max 
is higher for larger M. For given M, p max is higher for 
smaller J. A thermonuclear explosion is expected to oc- 
cur when the p max exceeds a certain critical value. Al- 
though the actual ignition density depends on the tem- 
perature (and thus M) , geometry and the explosion mod- 
els, we adopt 2 x 10 9 gcm~ 3 for th e ignition density as 
in the one-dimensional model W7 l|Nomoto et alJll984(l 
in the present approach with the complete degenerate 
model sequences. 

The criterion for the critical rotation of a rotating WD 
depends on the rotation law. In rigid rotation, which we 
adopt in the present study, it is given as 



where ij is the ratio of the equatorial centrifugal force to 
the gravity. A WD could break-up at equatorial plane 
when r\ becomes unity. 

Figure n show the area in the M — J plane where WDs 
can exist stably under these criteria. Each point repre- 
sents a WD of a different set of (p ma x, v)> which is trans- 
formed into (M,J). The sequence of WD models with 
pmax = 2-0 x 10 9 g cm~ 3 forms right edge of data area, 
and upper edge is composed of WDs whose ratio of sur- 
face gravity to centrifugal force are unity (77 = 1). As 
seen in Figure WDs of lower £ (i.e., with a smaller 
non-rotating region) have wider stable region in the pa- 
rameter space. For £ = 1, all WDs are on the M-axis 
(no-rotation spherical case). 

In Figures |21 and |3 examples of the structure of the 
accreting WDs just before the supernova explosion (i.e, 



p c = 2.0 x 10 9 g/cm 3 ) are shown. The equatorial radius 
and polar radius are normalized by the equatorial radius. 
In both figures, the left panels are almost nonrotating 
case (77 ~ 0, axis ratio q — 128/129; maximum value of 
q for rotating WDs in this method) and the right panels 
are almost break-up case (77 ~ 1). 

What would occur when a WD reaches the upper edge 
of the stable area by accretion? Is the WD broken up? 
Does it trigger the mass shedding and divides itself into 
sever al pieces? There is one possible answer for this ques- 
tion. [Paczvnski (1990) and iPopham fc Narav au (1990) 
found a solution that permits the WD to accrete without 
becoming unstable or triggering mass shedding, by com- 
puting a system consisting of a rigidly rotatin g WD and 
its vis cous accretion disk as one fluid, (see also lFuiimotol 

cnia)). 

Normaly, a rigidly rotating WD gains angular momen- 
tum from the Keplerian accretion disk due to viscosity. 
Inward angular momentum flow spins up the WD to in- 
crease ilo- When Slo exceeds the critical value (il cr ^) 
and rj exceeds unity, angular momentum is transported 
backwards from the WD to the disk while the mass of the 
WD increases. This state of rot ation is called "supercrit- 
ical rotation" ijPaczvnskill 1 990(1 . Deprived of its angular 
momentum, the WD slows down and when its angular 
velocity falls below O cr jf , the supercritical rotation ends 
and the WD gains both mass and angular momentum 
again. The viscous accretion disk acts like a stabilizer 
under the supercritical rotation. Then, if a WD reaches 
the critical rotation, we regard it will evolve along the up- 
per envelope in the J-M plane, and eventually explode 
as a SN la. 

2.3. Accretion Model 

In our accretion model, the initial state is defined 
only by the initial mass (Mj nt ). We consider 0.60M Q < 
Mi n t < I.O7M0 for the initial mass range of a WD. This 
is the range where C+O WDs in the close binaries can 
form l)IImeda et alJll999l e.g,). Therefore, the range of 
£ is 0.42(= 0.60/1.40) < £ < 1.0 in the slow angular 
momentum transfer. 

The accreting WD gains mass AM and angular mo- 
mentum AJ at each step. The value of AM is fixed 
throughout the sequence. The new mass M new and an- 
gular momentum J new are expressed as 

M now = M + AM (13) 
Jnew — J + A J = J + jAM (14) 
where j is specific angular momentum of accreting gas 
as given by eq. © • 

When the new mass and angular momentum are ob- 
tained, the nearest data point in the stable area is 
searched. Next, two dimension interpolation is made us- 
ing the surrounding four points in M — J plane because 
we cannot directly caluculate the WD that have the ex- 
act value of M n ew and Jnew The obtained new potential 
and the specific angular momentum j are used for the 
calculation of the next step of accretion. 

In the slow-transfer case, the value of £ decreases as a 
WD accretes matter and its whole mass increases. Be- 
cause we only have finite resolution in £ plane, we can- 
not calculate the model of exact value of £ in each step. 
Therefore, we use the data sets of WDs of the nearest 
£, and switch them as £ decreases. The area in M — J 
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0.2 0.4 0.6 0.3 1 12 1.1 \,f 0.2 0.4 0.6 0.8 1 12 11 1 .6 

M(M-.-. ill M(Msun) 

Fig. 1. — The total data area in M — J plane where WDs can exist stably under the assumptions mcnsioncd in the text. The left panel 
is for the case of true rigid rotation (£ = 0), and the right panel is for the case of slow angular momentum transfer (£ = 0.9), respectively. 




0.2 0.4 0.6 0.8 1 0.2 0.4 0.6 0.8 1 



Fig. 2. — Density contour in the meridional cross section for p c = 2.0 X 10 9 (just before the supernova explosion). Rotation law is true 
rigid rotation (£ = 0). The left panel is nearly spherical case (r) ~ 0) and the right panel is nearly critical rotation case (jj ~ 1). 




0.2 D.4 D.6 0.8 1 0.2 D.4 06 0.8 1 

Fig. 3.— Same as Figure 2, but for £ = 0.9. 
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plane where WDs can exist stably without the critical 
rotation depends on £, and WDs sometimes reach the 
critical rotation before £ becomes low enough to switch 
the data sets because of the low resolution in £ plane. 
We consider WDs evolve along the upper edge of each 
data sets same as the fast-transfer case until £ becomes 
low enough to use the next data set of £. 

3. RESULTS 

Many WD models are calculated for various parame- 
ters. For each initial mass, the evolution of WD is fol- 
lowed. 

3.1. Fast Angular Momentum Transfer Case 

Results for the fast angular momentum transfer case 
are shown in Figure As described in §2, £ = 
throughout the evolutionary sequence. The accreting 
WD reaches the upper edge of the data table regard- 
less of its initial mass as seen in Figure 0] It means that 
all WDs rotates fast enough to be affected by the su- 
percritical rotation. The mass when the accreting WD 
reaches the supercritical rotation (M er it) varies from 0.70 
to 1.21M0, as summarized in Table ETT1 Although the 
property of WDs is diverse when they reach supercriti- 
cal rotation, they will converge to a single point, which 
is at the upper right corner in the left panel of Figure 
and explode as supernovae. The final mass Mfi n is 
1.48Mq and the final total angular momentum Jfi„ is 
4.63 x 10 49 erg-sec. The effect of rotation contributes to 
the homogenity of SNe la rathar than causes the diver- 
sity if angular momentum transport is much faster than 
accreting timescale. 

TABLE 1 
INITIAL AND CRITICAL MASS 
FOR FAST TRANSFER 



M in t M crit 

(Ms) (Mp) 



J 

(10 49 erg-sec) 



0.60 
0.70 
0.80 
0.90 
1.00 
1.07 



0.701 
0.814 
0.925 
1.03 
1.14 
1.21 



5.80 
6.70 
7.44 
7.90 
8.08 
8.03 



In Figure |3] the radius and maximum density of the 
accreting WD is plotted against the mass as the WD ac- 
cretes matter from the companion star. Initial masses of 
the WDs in each sequence is 0.6, 0.8 and 1.07M Q , respec- 
tively. Before reaching the critical rotation stage, radius 
of the WD increases because the net gravitiy is decreased 
as the rotation grows faster. When supercritical rotation 
occurs and rotation is stabilized, the radius decreases as 
the WD accretes matter to increase its mass. For the 
central density, on the other hand, the sequence is not 
monotonic before WDs reach supercritical rotation. Un- 
der the supercritical rotation, the central densities in all 
sequences converge and increase toward SNe la. 

3.2. Slow Angular Momentum Transfer Case 

For the slow angular momentum transfer case, results 
are shown in Figure Because the value of £ changes 



as the WD accretes matter to increase its mass, the ef- 
fect of supercritical rotation suppresses the gain of angu- 
lar momentum continuously, although sequences are not 
smooth because of the limited resolution in parameter 
space of £. 

Unlike the fast angular momentum transfer case, the 
final state of WDs are not uniform but depend on their 
initial mass (Mnt)- The total mass and angular mo- 
mentum of WDs just before the supernova explosion 
(Mfin, Jgn) are summarized in Table J nn is smaller 
than the fast angular momentum transfer case, because 
the value of £ is non-zero and the region on the M — J 
plane where WDs can exist stably is narrower. As 
seen in Table 13.21 WDs with smaller M; n t have larger 
mass and angular momentum when their central density 
reaches 2.0 x 10 9 g/cm 3 . The difference of mass is small 
(~ 0.02M Q ); however, the difference of angular momen- 
tum is rather large (~ 2 x 10 49 erg • sec). 

TABLE 2 
INITIAL MASS AND FINAL STATE 
FOR SLOW TRANSFER 



M n t Mfi„ 



</fin 



£fin 



(M ) (M ) (10 49 erg-sec 



0.60 
0.70 
0.80 
0.90 
1.00 
1.07 



1.47 
1.47 
1.46 
1.46 
1.45 
1.45 



4.15 
3.94 
3.72 
3.43 
3.15 
2.76 



0.42 
0.49 
0.56 
0.63 
0.69 
0.75 



4. DISCUSSION 

We obtained the result that the accreting WDs in close 
binaries have the tendency that WDs with smaller ini- 
tial mass (Mi n t) have larger final mass (Man) and an- 
gular momentum (J nn ) when we assume that angular 
momentum tra nsfer is sufficientl y slow compared to the 
accretion rate ()Yoon et al]l2002j) . 

It is reported that brightnesses of SNe la show 
the following depe ndence on their parent galaxies 
<|Hamuv et al J 127)0(1 . The bright SNe la are observed 
only in spiral galaxies, while dimmer SNe la are found 
in both spiral and elliptical galaxies. Effects of rotation 
might qualitatively explain this tendency as follows. 

In elliptical galaxies, star formation has stopped long 
time ago, so the massive stars have already evolved. 
Therefore, a WD with M; n t that is too small cannot ac- 
crete enough mass to make a supernova explosion be- 
cause no companion star can supply enough matter to it 
(?). On the other hand, in spiral galaxies, star forma- 
tion is still active, so the companion stars of WDs can be 
massive. Therefore, WDs of small M- lnt can get enough 
mass to become SNe la. 

With the effect of rotation, accreting WDs with smaller 
Mi n t grow to larger M nn and J nn at the explosion. Al- 
though the difference of M nn is small, the difference of 
Jgn might affect the explosion and following nucleosyn- 
thesis. If there is a relation that 56 Ni is synthesized more 
for larger J, accreting WDs with smaller Mj n t grow to 
larger J nn and are expected to become brighter SNe la. 
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Fig. 4. — The evolutionary tracks of accreting WDs for fast transfer case (£ = throughout the sequence). The left panel shows the 
track of M- lnt = O.GMq and the right panel shows one of M- ln t = 1.07Mq. 





Fig. 5. — Radius and central density of accreting WDs are plotted against mass during accretion. The left panel shows is for radius and 
the right panel is for central density. M; nt is 0.6, 0.8, 1.07Mq. 





Fig. 6. — Same as Figure 4, but for slow transfer case. 
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This implies that the upper limit of SNe la brightness is 
higher in spirals than in elliptical galaxies. In this way, 
SNe la in spiral galaxies can originate from WDs with a 
wider range of M; nt , thus having a larger dispersion of 
brightness than those in elliptical galaxies. This could 
explain the observation that the most luminous SNe la 
appear to be observed only in spirals, while dimmer SNe 
la are observed in both spiral and elliptical galaxies. 

In the above discussion we assumed that angular mo- 
mentum transfer inside the accreting WDs are suffi- 
ciently slow compared to the accretion rate (M), and 
our calculation is hydrostatic. To examine whether these 
assumptions are valid, more precise study using multi- 
dimensional hydrodynamic caluculation is needed. As 
an accreting WD gains matter from their companion 
star, the whole binary system also evolves and M should 
change. Therefore we have to consider not only the ac- 
creting WD itself, but also the whole binary in order to 
treat angular momentum transfer precisely. 

Furthermore, we do not have precise theoretical knowl- 
edge of flame propagation and nucleosynthesis in the su- 
pernova explosion. In addition, the effect of rotation 
on these phenomenona is not clear. We have to exe- 
cute multi-dimensional explosion simulations to investi- 
gate whether the rotation causes the significant diversity 
of the luminosities of SNe la. 

5. CONCLUSIONS 

We have computed the structure of the rotating ac- 
creting WDs of various axis ratio (q) and central density 
(p c ). For the rotation law, rigid rotation is assumed and 
the ratio of the nonrotating core to the whole WD (£) is 
introduced. For each value of £, region on M — J plane 
where WDs can exist stably is obtained. When we sim- 



ulate the evolutionary track of WDs in M — J plane, we 
consider they continue to accrete matter u nder the effect 
of supercritical rotation ijPaczvn ski 1990) after reaching 
the critical rotation. 

The effect of rotation depends on the angular momen- 
tum transfer timescale compared to the accretion rate. 
If angular momentum transfer is sufficiently fast, all an- 
gular momentum brought to the WD is averaged and the 
WD rotates rigidly. WDs reach critical rotation and con- 
verge to the single state regardless of their initial state. 
Therefore, rotation contributes to the uniformity of SNe 
la when angular momentum transfer is fast. On the other 
hand, when angular momentum transfer is slow, there ex- 
ists a nonrotating central region of the WDs. The size of 
this nonrotating region depends on the initial mass of the 
accreting WDs, as the final state depends on the initial 
state. 

We have found that for accreting WDs, a smaller initial 
mass Mi n t results in a larger final mass Mg n and angular 
momentum J nn at explosion. This diversity of M nn is 
about 0.02M Q and may be too small to contribute to 
the observed diversity of brightness of SNe la. On the 
other hand, the difference in J nn could be sufficiently 
large to affect the speed of flame or nucleosynthesis in 
the explosoion. Further study is needed to demonstrate 
the effect of rotation on the explosion and whether the 
rotation can be the source of the diversity of brighnesses 
of SNe la. 

We would like to thank Sung-Chul Yoon and Nor- 
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